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A fully nonlinear, dynamically consistent 
numerical model for solid-body ship motion. 

I. Ship motion with fixed heading 

By Ray-qing Lin 1 ’* and Weijia Kuang 2 

1 David Taylor Model Basin, Carderock Division, NSWCCD, 

West Bethesda, MD, USA 

2 NASA Goddard Space Flight Center, Greenbelt, MD, USA 

In this paper, we describe the details of our numerical model for simulating ship solid- 
body motion in a given environment. In this model, the fully nonlinear dynamical 
equations governing the time-varying solid-body ship motion under the forces arising 
from ship-wave interactions are solved with given initial conditions. The net force 
and moment (torque) on the ship body are directly calculated via integration of the 
hydrodynamic pressure over the wetted surface and the buoyancy effect from the 
underwater volume of the actual ship hull with a hybrid fmite-difference/fmite-element 
method. Neither empirical nor free parametrization is introduced in this model, i.e. no 
a priori experimental data are needed for modelling. This model is benchmarked with 
many experiments of various ship hulls for heave, roll and pitch motion. In addition to 
the benchmark cases, numerical experiments are also carried out for strongly nonlinear 
ship motion with a fixed heading. These new cases demonstrate clearly the importance 
of nonlinearities in ship motion modelling. 

Keywords: ship motion; nonlinearity; hydrodynamics; solid body; numerical model 


1. Introduction 

For many decades, much effort has been devoted to modelling a ship’s motion 
at sea in naval architecture (Lewis 1989). Accurate prediction of ship motions 
in real time is critical for improving operability and preventing large-amplitude 
ship motions, e.g. vessel capsize. This is particularly important to today’s modern 
fast passenger ferries, high-powered naval vessels and cargo ships. 

The earliest effort on modelling ship motion can be traced back more 
than half a century. For example, St Denis & Pierson (1953) first proposed 
a method to predict the statistics of ship responses to a realistic seaway. 
Using spectral methods developed in applied mathematics, they established 
a relationship between the spectral density of ship responses and the input 
ocean wave spectrum. Later, Peters & Stoker (1967) developed a different 
algorithm for thin ship body motion: a hrst-order theory using a systematic 
perturbation procedure with the ship’s beam and unsteady motion assumed 
* Author for correspondence (ray.lin@navy.mil). 
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comparable (of similar magnitudes). This approach was further improved by 
Newman (1961) with a set of small parameters and more appropriate body 
boundary conditions. 

Later, the strip theory was developed by Ogilvie & Tuck (1969) in which linear 
ship motion coefficients are introduced with the slender-body assumption. These 
coefficients include added mass and damping coefficients used in heave and pitch 
motions. The surface integration used in this theory complicates its computational 
implementation. The theory is also inconsistent: the formulation is applicable 
in the short-wavelength domains, while the slender-body approximation works 
in the long-wavelength domains. An interpolation theory (Maruo 1970) and a 
unified theory (Newman 1978) were then proposed to reduce this inconsistency. 

Introduction of the more complex (but still linear) Neumann -Kelvin approach 
in numerical modelling greatly advanced the ship motion research (e.g. Dawson 
1977; Magee 1994; Shin et al. 1997), because of its capability to model arbitrary 
ship surfaces. However, it has its own limitations, e.g. in solving a forward-speed 
Green function in finite water depth (Beck & Reed 2000). 

When nonlinear terms are comparable to, or even stronger than the linear 
terms, new models are then needed. As an intermediate step to fully nonlinear 
models, several ‘blending methods’ were developed. For more details, we refer the 
reader to the ISSC report on extreme hull girder loading (ISSC 2000). 

A different approach was introduced by Wilson et al. (1998) and Gentaz et al. 
(1999) to solve the Reynolds averaged Navier-Stokes (RANS) equations in the 
time domain in which an iterative method is used for steady solutions, and a time- 
stepping method is used for unsteady solutions. The outcome of this approach 
is inconclusive, partly owing to insufficient numerical results. More seriously, it 
faces convergence problems for strongly nonlinear ship motion. 

Several other ship motion models were developed based on potential flow, 
e.g. Lin et al. (1986). However, they only partially include nonlinearities, and 
also use many empirical and free parameters, thus limiting their applications to 
ship design. 

In a more recent review by Beck & Reed (2000), approximately 80 per cent 
of ship motion models are based on the strip theory for its simple numerical 
implementation and its flexibility on ship hull forms. This implies also that most 
models are not suitable for strongly nonlinear ship motion. 

To model strongly nonlinear ship motion, it is therefore necessary to develop 
a methodology completely different from the above traditional approaches. In 
particular, it should be dynamically and mathematically consistent, numerically 
efficient, and independent of a priori experimental data. 

For this purpose, we have developed a new generation ship motion model. In the 
first phase of our effort (Lin & Kuang 2004, 2006; Lin et al. 2005), we developed 
a ‘steady ship motion model’ in which only the nonlinear interactions of surface 
waves (including ship motion-generated waves) are simulated. The response of 
the ship hull to these interactions (e.g. dynamic pressure force and the buoyancy 
force from the displaced water) is not included in that model. In the rest of this 
paper we call it the phase 1 model. 

To continue this effort, we expand the model to include an independent 
component for simulating six-degrees-of- freedom, solid-body motion under the 
forces calculated by the phase 1 model. A challenge in this effort is to accurately 
model these forces when the ship position (e.g. underwater volume) varies in 
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time; that is, the position variation is determined by the six-degrees-of-freedom 
ship solid-body motion, but the motion changes are determined by these forces. 
To address these technical difficulties, we need to develop a suite of new 
numerical methodologies for this component. In particular, these methods must 
be numerically convergent and efficient to model dynamic dependences of the ship 
position and ship motion. The combination of this component and the steady ship 
motion model is called the ‘Digital, Self-consistent, Ship Experimental Laboratory 
Ship Motion Model’ (DiSSEL). 

In this paper, we report for the first time the mathematical details and the 
numerical benchmark results of DiSSEL, focusing in particular on the solid- 
body motion component. The paper is organized as follows: the mathematical 
equations of ship motion are in §2. Section 3 provides the benchmark results 
with respect to experiments and independent numerical models. Discussions are 
given in §4. 


2. Mathematical model 

This section provides the details of the mathematical equations, numerical 
algorithms and formulation. The description is focused on the solid-body motion 
component. There is only a brief summary of the phase 1 model; its details are 
in Lin & Kuang (2004, 2006) and Lin et al. (2005). 


( a) Model reference frames 

To integrate the solid-body motion component with the phase 1 model, one 
must examine the dynamical equations in two different reference frames. One 
is moving horizontally with the mass centre of the ship, with the origin set at 
the mean free surface (hereafter called the model reference frame). The other is 
attached to the ship body, with the origin at the ship mass centre (hereafter called 
the ship reference frame). The first is used in the phase 1 model. The latter is 
convenient for calculating the ship solid-body motion. The moment of inertia of 
the ship and the numerical grid of the hull are invariant in this reference frame. 

Ship movement can be decoupled into a three-dimensional translational motion 
(u c + v c ) of the mass centre x c of the ship, and a three-dimensional rotation 
motion Q about x c . In this description, u c is the ship velocity vector moving 
in the calm water, and v c is the translational motion vector responding to 
surface waves, 

v c = v x x + v y y + v z z, 

where v x , v y . and v~ are the surge, the sway, and the heave motions, respectively. 
Similarly, the angular velocity vector £2 is 

£2 = Q x x + Q y y + Q z z. 

In naval engineering, £2 X , Q. y and Q z are called roll, pitch and yaw, respectively. 

In the model reference frame the i-axis is the ship heading direction (from 
the stern to the bow), the y-axis points to the port and the 2 -axis is upward. 
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The ship mass centre x c = (0, 0, z c ) moves vertically with the velocity 



Zl 


(2.1) 


i.e. the heave motion of the ship. In addition, the reference frame moves 
horizontally with the velocity u c + v^ = u c x + v x x + v y f. 

The ship reference frame is often convenient to solve for the ship solid body 
motions (e.g. Goldstein 1980). In this reference frame the position vector is 
denoted by 

x s = (xs, y s , Zs) T , 


where the subscript s implies the quantity defined in the ship reference frame, 
and the superscript T means the transpose of the vector. Therefore, the position 
vector x in the model reference frame and x s in the ship reference frame can be 
transformed from each other via 

x = A • x s + x c 

and x s = A -1 • (x — x c ) 



where A is the transformation matrix and can be expressed with the Euler 
angles 6 


6 = 8 x x + d y y + 8 z z where — = i2 




cos d z 

— sin 6 Z 

O' 


COS dy 

0 

sin 6 y “ 


'I 

0 

0 " 

and 

A = 

sin d z 

cos d z 

0 


0 

1 

0 


0 

cos d x 

— sin d x 



0 

0 

1 


sin d y 

0 

COS dy_ 


0 

sin d x 

cos d x 


(2.3) 

Other Euler angle definitions can also be used as long as no two consecutive angles 
are defined relative to the same axis (Goldstein 1980). In particular, A" 1 = A T , 
and its time variation is 


-A = flxA. (2.4) 

d t 


For a rigid ship its boundary xf does not vary in time. However, in the model 
reference frame, its time variation is given by 



d . , y , dA v dx c 
— ( A - if + !„) = — - if + — 


= (fi x A) • xf + v z z. 


(2.5) 


In DiSSEL, equations (2.1), (2.4) and (2.5) are used to update the ship 
boundary x“. 


( b ) Hydrodynamic equations of the surface waves 

The surface waves are solved in the model reference frame. The fluid is 
assumed inviscid and incompressible. Therefore, the fluid flow can be described 
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by a single velocity potential cp. and can be solved by the following nonlinear 
hydrodynamic equations: 


and 




dcp 

dt 


2 V< P“ (u c + v h ) 


_ V 

Vcp + g • z + - 

P 


d(u c + v h ) 
dt 


■ x — vV^cp = 0, 


(2.6) 


(2.7) 


in —H < z < ry (77 is the free surface elevation and H is water depth). In addition, 
there are three boundary conditions 



dy dcp 

T7 + (V hV) ■ [Vh<P - (u c + v h )] = — at z = 7 ), 

dt dz 

(2.8) 


dcp 

— = 0 at 2 = —H 
dz 

(2.9) 

and 

hx • V(p = hv • [u c + v c + £2 x (x s — x c )] for x e x“. 

(2.10) 


In equation (2.7), g is the gravitational acceleration vector, p the total pressure, V 
the gradient operator and V>, the horizontal gradient. In equation (2.10), riv is the 
normal unit vector of the ship surface X s . The acceleration term in equation (2.7) 

d(u c + v h ) 
dt ' X 

is from the acceleration of the model reference frame. 

The inviscid fluid approximation works well for large-scale motions. However, 
dissipation is important for small-scale flow (also called the sub-grid processes 
if their spatial scales are smaller than numerical resolutions), which arises from 
nonlinear interactions of large-scale waves, and/or from other turbulent processes, 
such as wave-breaking processes. To accommodate this, a small dissipative 
coefficient v is introduced in equation (2.7). It also acts to ensure numerical 
convergence with moderate resolution while retaining correct large-scale flow 
structures of interest. For more discussion on this dissipative effect, we refer the 
reader to Lin & Lin (2004). It should be pointed out that this approach has been 
commonly used in computational fluid dynamics for many decades (e.g. Spiegel & 
Veronis 1960; Yanai 1983). It was also successful in our earlier simulation (Lin 
et al. 2005; Lin & Kuang 2006). 

In addition to the boundary conditions (2.8), (2.9) and (2.10), we also need 
conditions on the numerical domain boundaries. In our model, the far-held 
radiation conditions (as x — » 00) are used on the borders of the numerical domain 
that is confined by the forward boundary Tf and the side and aft boundaries T a . 
Thus, the radiation boundary condition on Tf is 

^ = ^ + K + Vh), V = Ve (2.11) 

dx dx 

where cp e and rj e are the velocity potential and the free surface elevation of the 
environment (which vanish in calm water). The condition (2.11) implies that 
the waves generated by the ship do not radiate ahead of the ship (provided that 
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the ship speed is faster than the wave propagation speed) . Open boundary condi- 
tions are applied on S a that permit waves to travel freely across the boundaries 

V<p = V<p e + ( Uc + v h ), v = Th (2-12) 

where the over line represents the spatial average along the boundaries. 
Equation (2.12) also ensures conservation of mass for incompressible flow. 


(c) Dynamic equations for the ship solid-body motion 

A rigid ship motion can be decomposed into the translational motion v c of its 
mass centre x c and the rotation motion Q about x c (Goldstein 1980). The mass 
centre of the ship is defined by 


X c 


1 

m s 


p s xdE 


(2.13) 


where nig is the total mass of the ship, G s the total ship volume and p s the density 
of the ship. In particular 


dx c 
d t 


= u c + v c . 


(2.14) 


In the model reference frame, x c = (0, 0, z r ). and the variation of z c is given by 
equation (2.1). 

The kinematics of the rotation in the model reference frame is also very simple. 
Given any point x 1 on the ship boundary, we have 

^(x s -x c ) = £2 x (x s -x c ). (2.15) 

Equations (2.14) and (2.15) describe the six-degrees-of-freedom, solid-body 
motion of the ship. 

The dynamics of the ship motion is governed by the classical mechanics 
equations: 

to s ^ + T> v v c = F, (2.16) 

at 

and 

^(I.fl) + flx(I-O) + D 0 (I-fl) = r, (2.17) 


where I is the tensor of the moment of inertia of the ship; D v and Do are the 
integrated dissipative effects (e.g. drag) of the fluid on the ship hull, and F and 
r are the net force and the net torque (moment) on the ship hull, respectively. 
The net force (moment) includes contributions from the dynamic pressure on the 
ship hull, and from the displaced water 


and 


F = F p + F ff , 



npdT ■ , 


E® — (Pw Twet m s )g 


(2.18) 
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and 


r = r p + r\ 


r p = - 


dT(x s — x c ) x n p 


■^wct 


r p — (x we t x c ) X (iris p w V we t)g 


(2.19) 


In equations (2.18) and (2.19), p w is the water density, x wet the geometric centre 
of F we t, -Swet and V W et are the wetted surface and the underwater volume of the 
ship, respectively. The moment of inertia is 


lij — 


p s [(x - x,.) ;! rh, - (x - x c )j(x - x c )j ] d V, 


(2.20) 


(<5 is the Kronecker delta function). 

It should be pointed out that the two dissipative coefficients, D v and Dq, are 
from various sources, such as friction (ignored in the large-scale potential flow 
approximation) owing to small-scale flow generated by, for example, bilge keels 
and wave-breaking. In many cases, the dissipation owing to the bilge keels is 
very strong, sometimes reaching approximately 20-25% of the net moment on 
the ship. The small-scale flow, and thus the dissipation, depends on the incident 
waves and the ship profiles as well (Lin & Kuang 2008). Friction from wave- 
breaking is in general weak and contributes only approximately 1 per cent of the 
ship acceleration (Lin & Kuang 2007). 

In numerical modelling of the ship motion accurate computation of the forces 
and torques (moments) is critical. This depends on several important factors. 
By equations (2.18) and (2.19), we can observe that accurate determination of 
the forces and the torques (moments) depends on two issues: an accurate wetted 
surface T we t (and therefore the underwater volume and an appropriate 

numerical algorithm for surface integration. In addition, the dissipative effect 
must also be modelled consistently. 

The wetted surface T we t varies in time, and can be determined if x c and 6 are 
known. They can be determined if F and T are given, but F and T depend on T wet 
via equations (2.18) and (2.19). Therefore, the relationships among the system 
variables are very nonlinear. In other words, all degrees of ship motion couple to 
each other. 


(d) Numerical methods 
In DiSSEL, the numerical state vector S is 

S = (<p,'n,v c ,Q,6,z c ). 

A hybrid scheme is used in the model to solve S. In this method, the velocity 
potential q> and the free surface elevation p are approximated by finite Fourier 
expansions, e.g. 

<p(x, y, z, t) = ^ <p(z, t)e l( ' kmX+kaV ^ + c.c. 

n,m 

The ship surface xf is described by either a finite-element mesh (unstructured 
grid) or a finite-difference mesh (structure body fix grid), and is fixed in the ship 
reference frame. 
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The numerical workflow is as follows: First given the initial state vector S(<o), 
the nonlinear wave-wave interactions are solved via fast Fourier transform (FFT) 
transform between the spectral space and the physical space. The wetted surface 
T we t at to is determined from S (to). Then the forces F and the torques T are 
evaluated via equations (2.18) and (2.19). After those calculations, a third-order 
Runge-Kutta scheme is used to update the state vector S (to + At). 

The advantage of the hybrid method is obvious. In the ship reference frame, 
the mesh is fixed and invariant, and is therefore determined before the time 
integration. Updating of (p and tj is carried out by the pseudo-spectral algorithm. 
Therefore, the central processing unit (CPU) time (flops) at each time step is 
determined by the FFT transforms, which is of the order (A log A) 2 . Compared 
with that of A 4 for the regular finite-difference scheme, this is computationally 
much more efficient. Another advantage of the hybrid method is its numerical 
accuracy. The ship hull is much smaller than that of the entire numerical domain. 
The grid sizes for surface wave calculations are not appropriate for the numerical 
integration of equations (2.18) and (2.19). However, the finite element /difference 
mesh x/' can provide sufficient resolution for the surface integration. 

Update of x E is determined by equation (2.5), i.e. by z c , £2 (and therefore A). 
Together with the updated free surface elevation 17, we can update the wetted 
surface T wet . Update of the normal vector n 2 can be calculated either from the 
updated x s , or more efficiently, directly from 


dri\ 

d t 


= £2 x n £ . 


Convergence of the numerical solutions is discussed in detail in Lin et al. 
(2005) and in Lin & Kuang (2006). Details on the integrated dissipative effects 
in equations (2.16) and (2.17) are described in Lin & Kuang (2007, 2008). 


3. Model benchmark results 

To validate the DiSSEL, we have carried out simulations for many cases with 
either experimental data or results from independent numerical ship motion 
models, e.g. LAMP (Large Amplitude Ship Motion Programme; Lin et al. 1986). 
In these cases, the ship motion generated wave patterns, nonlinear wave profiles 
around the ship, pressure distributions and unsteady ship motions (roll, pitch 
and heave) are compared and cross-examined. Reported here are selected cases, 
which are the extreme scenarios of the benchmarked examples. A fixed heading 
is assumed in all benchmark cases. Reported cases include ship motions in 
Sea State 3 (SWH = 0.9144 m, Tp = 10.35s), Sea State 4 (SWH= 1.9m, Tp = 
10.75s) and Sea State 5 (SWH = 3.25m, Tp= 11.85s; note: SWH, significant 
wave height; Tp, wave period). 

(a) X -craft catamaran 

The X-craft catamaran is an example of a high-speed ship. Its Froude number 
reached Fr= 0.77 in experiments, and thus is appropriate for testing the model’s 
capabilities for very fast ship motion (one of the challenges for ship motion 
modelling). In this benchmark case, we simulated ship motion with three different 
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Figure 1. Pitch motion (O y ) of X-Craft in Sea State 4 head seas and Pr= 0.77. Solid line, DiSSEL; 
dashed line, data; circles with dashed line, LAMP. 



Figure 2. Numerical X-Craft ship wake in calm water by DiSSEL, Fr is Froude number. 
(a) Fr= 0.23, ( b ) Fr= 0.38, and (c) Fr= 0.77. 


Froude numbers: Fr= 0.23, 0.38, 0.77. Our results are then compared with 
experimental data and simulation results from an independent ship motion 
model LAMP. 

Figure 1 shows the results of the pitch motion in head seas for Sea State 4 and 
Fr= 0.77 from DiSSEL (solid line), experiments (dashed line) and LAMP (circles 
with dashed line). Compared with the LAMP results, our simulations are much 
closer to the experimental data. In particular, the amplitudes of the two sets of 
results are comparable. Better agreement is found for smaller Froude numbers 
as does the agreement between LAMP and experiment improve in these cases. 
Similar conclusions apply to the roll motion. 

We also plot in figure 2, the ship motion-generated wave distributions for the 
three Froude numbers. Though there is no experimental measurement for the 
cases, there are results from other ship motion models (only for small Froude 
numbers), which our results agree well with (Lin et al. 2005; Lin & Kuang 2006). 
The wave distributions in figure 2 demonstrate clearly the robustness of our ship 
motion model for high-speed ships. 
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Figure 3. Pitch motion (6 y ) of HSS at Sea State 5 heading 120° and Fr= 0.51. Solid line, DiSSEL; 
dashed line, data; circles with dashed line, LAMP. 



Figure 4. Similar to figure 3, but for the roll motion ( 0 X ). Solid line, DiSSEL; dashed line, data; 

circles with dashed line, LAMP. 


(b) HSS trimaran 

The HSS Trimaran is another type of high-speed craft. We simulate motion 
of this craft in bow seas (120° between ship heading and incident waves) with 
Sea State 5 and Fr= 0.51. In this case, there are both roll and pitch motions 
and it is thus ideal for understanding (nonlinear) interactions between different 
components of ship solid-body motion. 

The results for the pitch motion and the roll motion are shown in figures 3 
and 4, respectively. The line types for different results are defined similarly to 
those in figure 1. 

The simulation results (from DiSSEL and LAMP) of the pitch motion are 
similar and agree well with the experimental data, as shown in figure 3, but 
the roll motion results are different. The time series from DiSSEL agree well 
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Figure 5. Heave motion (tj z ) of PSD at the planning regime (Fr= 0.536) in Sea State 3 of irregular 
head sea. Solid line, DiSSEL simulation; dashed line, experimental data. 


with that of the experiment (in both amplitude and frequency). However, the 
agreement between that of the LAMP and the experiment is poor and becomes 
satisfactory only for smaller Froude numbers ( Ft < 0.4). The benchmark results 
indicate that accurate simulation of roll motion is much more difficult than that of 
pitch motion. This is because the roll motion depends strongly on nonlinear ship- 
wave interactions. These interactions become stronger as the forward speed of the 
ship increases. The good agreement between DiSSEL and experimental results for 
both pitch and roll motions demonstrates that these nonlinear interactions are 
properly implemented into DiSSEL. 


(c) Planning craft, propulsion system demonstrator 

Ship motions of the propulsion system demonstrator (PSD) in the high- 
speed planning regime are an example of a strongly nonlinear ship motion 
problem. In particular, the craft’s underwater volume is reduced substantially 
from the displacement regime (low speed) to the planning regime (high speed). 
Consequently, maintaining dynamical balances accurately at different stages is a 
numerical challenge. Since the balances determine the trim angle and the vertical 
displacement, numerical errors in these balances will lead to erroneous and often 
divergent solutions. Because of the difficulty, there has not been any ship motion 
model capable of modelling this craft in the planning regime. We select this 
benchmarking case to demonstrate the capabilities of DiSSEL to simulate strongly 
nonlinear ship motion. The mean values of the vertical displacement and trim 
angle related to the ship forward speed are approximately equal to those in calm 
water (e.g. Lin & Hoyt 2007), and the heave and pitch motions in the following 
figures only include those arising from the ship-wave and wave-wave interactions. 

Figures 5 and 6 show the heave and pitch motions at Fr= 0.526 (the planning 
regime) and Sea State 3 in head seas, respectively. The numerical time series (the 
solid line) agrees nearly perfectly with the experimental data (the dashed line) 
for both pitch and heave motions. 
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Figure 6. Similar to figure 5, but for the pitch motion ( 8 y ). Solid line, DiSSEL simulation; 

dashed line, experimental data. 



Figure 7. Similar to figure 5, but with Fr= 1.56. Blue line, simulated heave motion; pink line, 

incident wave height. 


In addition, we have carried out simulations with even higher designed planning 
speeds (Fr< 1 .56) . The results for Fr= 1.56 and Sea State 3 are shown in 
figures 7 and 8. There are no experimental measurements for benchmarking at 
these conditions. However, the numerical results are reasonable because the heave 
motion follows the incident waves with a small phase lag (as shown in figure 7), 
consistent with the empirical conclusions from experiments. Therefore, the results 
here can be used as ‘prediction’ for future experiments. 

In figures 9 and 10, we plot the force/moment (from the dynamic pressure) and 
the restoring force/moment (from the buoyancy) from the simulations. From these 
figures, we can observe that the two kinds of forces (moments) are comparable in 
magnitude, but nearly opposite in direction, indicating the establishment of the 
dynamical balance in the numerical solutions. 
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Figure 8. Similar to figure 6, but with Fr= 1.56. Solid line, simulated pitch motion. 



time (s) 

Figure 9. The numerical heave force F? and the restoring force F® (in the unit of 10 8 N) for the 
benchmark case in figure 7. Blue line, Ff; pink line, Ff. 


Planning craft sometimes exhibit low freeboard and bow submergence, often 
called ‘Plow-In’, as shown in figure 11 by the DiSSEL simulation. Plow-In occurs 
when the frequency of the ship motion generated waves is similar to the ship 
natural frequency. Numerical simulations show that Plow-In occurs when the 
forward speed is between 2.83ms -1 (Fr— 0.27) and 3.24ms -1 (Fr=0.31). This 
condition is very similar to that observed in experiments (Hoyt & Lin 2007). 
We want to emphasize here that this is the first successful numerical simulation 
of Plow- In. 
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Figure 10. The numerical pitch moment F^ and the restoring moment Fy (in the unit of 10 8 Nm) 
for the benchmark case in figure 8. Blue line, Fy\ pink line, Fy. 



Figure 11. Simulation of the Plow-In phenomena of PSD with U c = 3.1 m ( Fr= 0.2973) by DiSSEL. 


( d ) ONR hull form 

Ship motion is also sensitive to the hull geometry, under and above the water 
(defined in calm water) . The motions of two ship hulls with identical underwater 
geometry (in calm water) can be very different if the above-water geometries are 
not the same. This difference is significantly amplified in the resonant state in 
which the incident wave frequency is the same as the ship natural frequency, as 
shown in experimental results (figure 12). Simulation of the ship motion therefore 
requires correct evaluation of the forces and the moments on motion-dependent 
(thus time-varying) wetted ship surface. 
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A 

Figure 12. Normalized roll motion ( 6 x /ka ) of the ONR series Tumblehome hull and Flare hull for 
different wavelengths X (scaled by the ship length) in beam seas. Solid line, DiSSEL for Tumblehome 
hull; dashed line, DiSSEL for Flare hull; squares, data for Tumblehome hull; circles, data for 
Flare hull. 


To demonstrate DiSSEL capabilities of resolving forces and moments with 
time- varying wetted surfaces, we simulate the roll motions of the ONR series 
Tumblehome hull and Flare hull. Both have the same underwater geometry (in 
calm water), but the above water geometries are very different; one is slanted 
inward and the other is outward. Figure 12 shows the results of the normalized 
roll motions (scaled by the wave steepness) of the two hulls. From the figure, we 
can observe that the numerical solutions agree very well with the experimental 
data. We can also observe that the difference in the normalized roll motion of the 
two hulls reaches a maximum at resonance. Obviously, one could not obtain the 
correct results if the underwater geometry is assumed invariant (as often used in 
the linear ship motion model). 


4. Discussion 

In this paper, we provide an overview of the DiSSEL ship motion model, and 
the detailed description of the ship solid-body motion component. Included in 
this model are all wave-wave and ship-wave interactions. Consequently, all six 
components of the ship solid-body motion are dependent on each other and 
any change in one component will affect the rest via nonlinear interactions of 
the system. In addition, dissipative processes from the blocking effect (Lin & 
Kuang 2008) and from sub-grid effects (e.g. turbulence, wave-breaking, subscale 
processes) on the ship motion are included in the dynamical equations (2.16) and 
(2.17). In particular, the dissipative effect is formulated based on fundamental 
physics and mathematical properties. There is no other parametrization (e.g. 
added mass) or linearization applied in our model. With our approach, a single 
velocity potential is needed to model all nonlinear interactions. 
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The numerical algorithm used in our model is a hybrid of pseudo-spectral and 
hnite-difference/fmite-element algorithms, capable of resolving very complex ship 
hull geometries and broad surface wave spectra. The algorithm is computationally 
efficient and robust (Lin & Kuang 2006). 

Many benchmarking tests have been performed with the DiSSEL model. 
Reported in this paper is only a small fraction of them, which are chosen for 
the following reasons: (i) well-documented experimental and/or independent 
numerical simulation results; (ii) strongly nonlinear effects in ship motion; 
and (iii) numerically challenging problems. We also provide new simulation 
results (case c with Fr= 1.56) for future benchmark and prediction purposes. 
In particular, we would like to point out that the DiSSEL model can be used 
to simulate ship motion with any arbitrary forward speed, i.e. with the Froude 
number values far larger than those reported in this paper. 

In all benchmark cases our simulation results agree well with experimental 
data, while the agreement of independent models varies. In general, our simulation 
results outperform those of independent models. 

Our work is a step further towards providing a realistic ship manoeuvring 
model. By integrating the ship solid-body motion component and the surface wave 
model, we are able to determine ship motion, the motion generated waves and the 
responses to ship-environment interactions concurrently. However, the current 
version of DiSSEL is not complete. Future development includes simulation of 
ship motion for a ship manoeuvring in a seaway. 
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